THE PROBABILITY OF RECTANGULAR UNIMODULAR 
MATRICES OVER ¥g[x] 

XIANGQIAN GUO AND GUANGYU YANG 

Abstract. In this note, wc compute the probability that a k x n matrix 
can be extended to an n X n invertible matrix over Fq[x], which turns out 
to be (1 - q''-"){l - q''-^-") ■ • • (1 - g^-"). Connections with Dirichlet's 
density theorem on the co-prime integers and its various generalizations are 
C 2 ^ ■ also presented. 
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1. INTRODUCTION AND MAIN RESULTS 

"as 



For any set A of positive integers, its natural density is defined 



P< : (1) DiA) := litrr '^^ ^ ^^' ^l' ' ' ' ^>' . 



provided the limit exists, where | • | denotes the cardinahty of the corresponding set. 
Dirichlet [2] discovered an interesting density theorem that asserts the probabiUty 



2 ' of two integers to be co-prime is G/tt^, that is, 

►-->. . where gcd(TO, n) denotes the greatest common divisor of m and n, and C,{s) is the 

[^ ' Riemann's zeta function. Moreover, the probabihty of n integers to be co-prime is 

^. \{{mi,...,m„) g N" : 1 ^ mi,...,mn s^ iV,gcd(mi, . . . ,mn) = 1}| _ ^^ 

^^ I Kubota and Sugita [5] gave a rigorous probabiUstic interpretation to Dirichlet's 

density theorem. For the deep links between probability theory and number theory, 
please refer to Tenenbaum [12], Kubilius [4] and Kac [3]. 

Denote by Mkxn{R) the set of all fc x n matrices over a ring R. A matrix 
K^ \ A G Mkxn{R) is called unimodular if it can be extended to an n x n invertible 

'V^ ■ matrix. Using the concept of unimodular, we can give a re-statement of Dirichlet's 

C^ I density theorem and its generalization: the probability that a random 1 x n integer 

matrix is unimodular is C,^^(rL). Naturally, we can consider the matrix form of the 
Dirichlet's density theorem: the probability that a random fc x n (1 ^ fc ^ ri) integer 
matrix is unimodular? The key is to define the analogous natural density for the 
sets of the integer matrices. Maze, Rosenthal and Wagner [7] studied this problem 
and obtained that the probability that a random fc x n matrix is unimodular is 
n'l^n-fc+i C~^(j)- Furthermore, [7| also gave some interesting historical remarks 
for Dirichlet's density theorem. 

Then it is natural to consider the same question for matrices over polynomial 
rings of fields. As remarked in [7], the concept of natural density does not extend 



2000 Mathematics Subject Classification. 11C99, 15B33, 60B15. 

Key words and phrases. Natural density, unimodular matrices, Riemann's zeta function, q-zeta 
function, finite field, polynomial ring. 

1 



2 X. GUO AND G. YANG 

naturally to the ring F[a;] for a general field F. Fortunately, this can be done for 
the finite fields. 

Let Fq be a finite field consisting of q elements, and Vq[x\ be the polynomial ring 
over ¥q. To enumerate Fg[a:], let S be the set of all vectors a = (oq, oi, • • • ) with 
aiG{0,l,---,g — f} and a^ = for sufficiently large i. Then there is a one-to-one 
map X : S ^ ^+ = N U {0} defined by x(«0; o-ir ' ') = Si^o ^i^^ ■ For convenience, 
we denote a„i ~ X^^im) and TOq, = xipi)- Then for all m G Z+, we set 

oo 
fra{x) := ^ fljX', with Um = (oo, Ci, ■ • • ). 
1=0 

By a rigorous probabilistic method, Sugita and Takanobu |10| determined the 
probability of two polynomials over F^ to be co-prime when q is a prime, that is, 

,.. ,. |{(m, n) g {0, 1, . ■ . , iV - 1}2 : gcd(/^, /„) = f }| 1 

(4) hm —: = f . 

Af^oo N^ q 

Recently, several authors considered the question of the probability of n polynomials 
to be co-prime in Fg [x] . Morrison [5] introduced a concept of natural density 

(5) 5(5):=lim^^, 5 C Fjx], 

where Fd is the set of all the polynomials in F[a;] with degree no more than d. He 
first treated this problem in a heuristic way. Then through an elegant as well as 
rigorous argument, he deduced the following formula, 

(6) 

|{(/i,...,/„)6F^[x]:gcd(/i,...,/„) = l,deg(/,)«;d,f ^*<n}| ^ 1 

hm — rrz „ , , ; ttt rr; — 1 



|{/GFg[x]: deg(/)^d}|" 9"-i- 

Benjamin and Bennett [T] also studied this problem; their methods are much inter- 
esting using the Euclidean algorithm which is one of the oldest algorithms. Essen- 
tially, they also got the same conclusion. 

Our object in this paper is to extend the main results in [7] to the polynomial 
ring ¥q [x\ for any prime power g, that is, the matrix form of the Morrison's theorem 
dSl). Denote F = F, and ¥[x\ = ¥q[x] for short. Let M := A4x,i(F[x]) be the set of 
all A: X n matrices over ¥[x] and AIat be the subset of Ai consisting of all matrices 
with entries in {/o(a;), /i(x), ■ • • ,/Ar(a;)}. For any subset S* C A^, we define the 
natural density of S in A^ to be 

(7) DiS) := lim ^i^, 

^ ' A^oo \Mn\ 

if the above limit ([7]) exists. Note that, in the special case k — 1, the limit in ^ is 
just the limit of a subsequence in ([7]) with N ~ q'^. Thus our definition of natural 
density is a more rigorous one. 

To describe our result, i.e., to determine the natural density of fc x n unimodular 
matrices over F[x], we introduce the following g-zeta function 

oo 

(8) ^^(^■) ■■= 11(1 - -j^)-^ = n (1 - ^)-""' 

/ "^ m = l ^ 

where / goes through all irreducible polynomials in ¥[x\ and Lpm is the number of 
irreducible polynomials in F[a:;] with degree m. In the above statements and what 
follows, by irreducible polynomials we always mean monic irreducible 
polynomials (the constant polynomial f is excluded as usual). Now we are 

in the position to give our main result. 



THE PROBABILITY OF RECTANGULAR, UNIMODULAR MATRICES OVER F,[2;] 3 

Theorem 1. Let E be the set of all k x n unimodular matrices over F[a;], then the 
natural density of E is 

(9) DiE) = ]jii--—^)^i[ii^—r-^ n Cw- 

/ ^ m=l ^ j=n~k+l 

Remark 2. In other words, the probability that a k x n polynomial matrix can 
be completed by an {n ~ k) x n matrix into an n x n invertible matrix over ¥[x] 
is ri?=n-fc+i Cn'^'^O)- And via introducing the q-zeta function Cg(s); "^^ obtain a 
similar formula to the one in [7] for integer matrices 

j^n—k-\-l 

where E denotes the set of all k x n unimodular matrices over Z. 
Remark 3. There is an interesting equation 

oo ^ 

(10) II(i-^'°^*'^r' = E'^'^' = T^' 

where f is over all irreducible polynomials. For an interpretation of this equation, 
see [8| . Putting t = q~^ , j G N m ilO\) . we get 

(11) C-(^-) = l-^- 

Hence, we obtain an explicit formula for the natural density of E, 

(12) D{E) = (1 - -^){1 -^) • • • (1 - ^), 
and an interesting identity concerning {ipm}, 

oo 

(13) 11(1-^)^'" -1-^' V,EN. 

Remark 4. In particular, if we put k = 1, Theorem{l\ shows that the probability 
that n polynomials over the finite field ¥ to he co-prime is 



oo ^ 

C\n)= 17(1 



\vm ^i L 



n7n ' ^n— 1 

rn=\ ■=■ ^ 

This extends Dirichlet's density theorem (0j and yields the same conclusion as 
Morrison [5] and Benjam,in and Bennett [I] . However, noticing that we have defined 
the concept of natural density in a different way from them, as we mentioned before, 
one may find that our result actually is stronger than theirs. 

We end this section with some remarks on the square case k = n. One may 
notice that our proof of Theorem [T] in [2] docs not apply to this case (also, Cg(l) is 
not well defined in this case). However, we have the following result. 

Proposition 5. The natural density of n x n unimodular matrices over¥q[x] is 0. 

Remark 6. The proof of Proposition\^is quite simple and similar to the proof of 
Lemma 5 in [7], so we omit the details. If we make the convention that 0^(1) = 0, 
Theorem\^ also holds for k = n. 
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2. Proof of Main result 



In this section, we shall give the proof of Theorem [T] Before this, we need some 
preparations. 

Let P be a finite set of irreducible polynomials in F[a;] and V be the set of all 
irreducible polynomials in F[a;]. Denote by E-p the set of all matrices A € Ai = 
MkxnC^lx]) such that the gcd of all full rank minors (i.e., minors of rank k) is 
co-prime to all elements in V. Recall that E is the set of all unimodular matrices 
in Ad, that is, each matrix in E can be extended to an n x n invertible matrix. It is 
well known that A G A4 is unimodular if and only if the gcd of all full rank minors 
is 1 (See [3], [H] and [13] for more details). Thus we have E = f]^ E-p. 

Lemma 7. Let E-p be defined as above, then we have 

fe-i ?i 

(14) D{Ep)^Y[Y[{l-q^'-"')'^''sU)^^ Yl Y[{l-q-''^°^^f^). 

j=o fev j=n~k+i fer 

Proof. Denote fp := Yifevf ^^"-^ '^'P ~ cleg(/-p). For any positive integer iV, 
consider the maps, 

(15) t:: Mn^ A4x„(F[a;]/ (fp)) 

and 
(16) 

i; : A/fcx„(FM/ {fp)) -^ Ah.n ( J] I^N/(./) ^ H ^'^''x« i^^^/if)) ' 

\f£V J f£V 

where {fp) and (/) denote the ideals generated by fp and / respectively. Note 
that TT is the canonical projection via modulo fp and tA is a composition of two 
canonical isomorphisms as F vector spaces. 

Taken any A G j\4n- It is easy to see that A G Ep if and only if the com- 
ponent oi ip o '^{A) in Mkxn{V[x]/{f)) is full rank for each f G V. Since f G V 
is irreducible, we know that F[a;]/(/) ~ F^dGgt/). Let F^daeif) denote the set of 
all full rank k x n matrices over the finite field F^dog(/) . It is well known that 

l-Fgdogc/) I = ]TjZo{<l"'^°^^^^ - q^'^^s^-'^') (See Ej page 455, for example). 

First suppose that N = mq'^'^ — 1 for some m G N. Then it is easy to see 

{fi{x) : O^l^N}^ {fs{x)x''^ + ft{x) : O^s^m- 1,0 ^t^q""^ -1}. 

For any fixed Sj s ^ ?7i — 1, the following projection is one-to-one: 

{fs{x)x''^ + ft{x) ■.O^t^q'i^-1}^ mx]/{fp) 

and the canonical projection 

{fi{x): 0^l^N}^¥[x]/{fp) 

is TTi-to-one. As a result, the projection map tt in (jlSp is 77i*''"-to-one. Thus we 
obtain that 

\EpnMN\--m'''' ■ 11 |Fgdcg(/,t 
fev 
fc-i 

= ™''" ■ n n q"'^°'^^^H^ - g(^'-")'^°s(/)) 

fev 3=0 

fc-i 
(17) = {mq'^^)''" ■ Yl ]J(1 -q(j-")'i'=s(/)) 

f&V3=0 
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Now suppose that N is any positive integer. There exists some m,r E Z+ such 
that iV + 1 = mq'^'P + r with ^ r < g'*''. For convenience, set TV := viq'^'^ — 1. 
Then by the definition of the natural density, we have 

D[E^)= hnr l^-^^^' 



N^oo \M 



N\ 



(18) = hm 



\Ev f^M^\ + \Evf^ {Mn - M^) 



Note that \Mn - Mf^\ < rkn{N + 1)''""\ which gives that 



\EvniMN-M j^)\ 



(19) lim ^^T^^^^TTi;^^ = 0. 



Thus, from dUl), dUl) and ^, we obtain 



D{E-p) = hm 



\EvnMff\ 



N^oo {N + 1)'^" 

AT^oo {N + I)''" 

feVj=o 

n 
j=n-k+l feV 

This completes the proof of Lemma [T] D 



Proof of Theorem [TJ Suppose that k < n. For any irreducible polynomial 
/ G Fg [x] , let Hf C A^ be the set of all matrices whose gcd of all full rank minors 
is divisible by /. Denote qf :~ q'^^&if) ^ then by Lemma [71 

D{Hf) = l-D{Ey}) 

<i- 1- E ly 

y j=n-k+l 

j=n-k+l 9/ W~l) 9/ 

Let Vt be the set of all irreducible polynomials with degree no more than t and 
denote Et ~ E-p^ . Since 

{Et\E)C y Hf, 
fev\Vt 
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we have 



\iEt\E)nMN\ ^,. \i[Jfev\v,Hf)nMN\ 
lim sup — — — j ^ lim sup -—— — ■ 

N^oo \Mn\ N^oo \Mn\ 

J2fev\Vt \HfnMN\ 



^ lim sup ■ 

^ y lim sup 



HfCiMNJ 



Sev\Vt fevyPt ^ 



Z^ «2dcg(/) Z^ 



2(p„ 



«2dcg(/) Z^ q2m ■ 

/e-p\-Pt m=t+l 

It is well known that all irreducible polynomials with degree m can divide x' — x, 
which has no multiple roots. Thus mtpm ^ 9™ and 



\{Et\E)nMN\ ^ y 2 



hm sup 7— — i s$ > 

AT^oo \Mn\ J—'^, mq 



< E ' 



Now note that E n Mn '^EtH Mn and £■ n Mn ^ EtH Mn - (Et \E)n Mn, 
which imply that 

^. .JEnMNl^y . r\Etr\MN\ ,. \iEt\E)nMN\ 

lim ml — — — — ■ — ^ hm mt — — — — ■ lim sup ■ 



A^^^^ \Mn\ N^^ \Mn\ JV->oo \Mn\ 

g*(g~l) 



and 



lim sup — r-rz — : — ^ hm sup — r— - — ; lim mi r— - — ; 

A^oo \Mn\ A^oo \Mn\ N^ca \Mn\ 

^ D{Et), 

for all t € N. Let t tend to oo, and recall that the number of monic irreducible 
polynomials with degree m in F[a;] is ipm- From Lemma [7] and the definition of 
g-zcta function, we can conclude that 

lim l^,^;^,"' = lim D{E,) 

A^oo \Mn\ *^°° 



lim n TT (l-q-J'*°e(/)) 



t—^oo 






hm n n (1 

j—n — k^l m—1 
n +00 n 

= n 11(1-5""")'^"= n CO')- 

j—n — k-\-l 711 — 1 j—n — k-\-l 

This completes the proof of Theorem [TJ D 
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